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In the present work the primordial Big-Bang Nucleosynthesis (BBN) and weakly interacting mas- 
sive particle (WIMP) dark matter are discussed in a certain class of modified gravitational theories, 
namely f(R) ~ R" gravity. The new gravitational model is characterized by a single parameter n. 
First we determine the conditions under which the theoretical predictions for the 4 He abundance are 
in agreement with the observations. More precisely, during BBN the physics is known and all the 
parameters are known. The only free parameter to be constrained is the power n related to the new 
gravitational model. After that, for cold dark matter we use the value of n determined from the BBN 
considerations and determine how the mass of the dark matter particle is related to the annihilation 
cross section in order for the cold dark matter constraint to be satisfied. 
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I. INTRODUCTION 

There is accumulated evidence both from astrophysics and cosmology that about 1/4 of the energy 
budget of the universe consists of so called dark matter, namely a component which is non-relativistic 
and does not feel the electromagnetic nor the strong interaction. For a review on dark matter see e.g. QjJ] - 
Although the list of possible dark matter candidates is long (for a nice list see e.g. (2J]), it is fair to say that the 
most popular dark matter particle is the LSP in supersymmetric models with R-parity conservation [3]. The 
superpartners that have the right properties for playing the role of cold dark matter in the universe are the 
axino, the gravitino and the lightest neutralino. By far the most discussed case in the literature is the case of 
the neutralino (see the classical review probably because of the prospects of possible detection. On the 
other hand, primordial Big-Bang nucleosynthesis (BBN) is one of the cornerstones of modern cosmology. 
In the old days, BBN together with Hubble's law and CMB supported and strengthened the Hot Big-Bang 
idea. Nowadays, BBN can be used to test and constrain possible new physics beyond the standard model. 
The new physics may be either due to exotic particles predicted by particle physics model or due to a new 
expansion law for the universe predicted by a new gravitational model. For a recent review on BBN see 
e.g. [5]. In the present work we shall be interested in a class of new gravitational models of the form 
f(R) ~ R n , where the power n is the only parameter that characterizes this class of models. Although in the 
literature the authors usually discuss this kind of modified gravitational models in the late times universe 
(see e.g. [(], 7]), here we wish to discuss this class of models in the early universe. In iQ] the authors were 
interested in the baryon asymmetry in the framework of gravitational baryogenesis proposed a few years 
agolB. 

In this Letter we wish to study this class of gravity models in two respects, namely primordial Big-Bang 
nucleosynthesis (BBN) and WIMP dark matter. Our investigation will allow us to first derive the allowed 
range for the power n, and to see how different this class of models can be compared to general relativity. 
Then for these values of n we determine how the WIMP mass has to be related to its annihilation cross 
section so that the cold dark matter constraint is satisfied. As a matter of fact, already in QSQ, the authors 
have mentioned that obtaining the right baryon asymmetry in agreement with BBN requires a value of n 
close to unity. Their discussion was based on the argument that the temperature relevant for BBN should 
be within the range 0.1 — 100 MeV. Here, however, we perform a more accurate investigation by actually 
computing the cosmological helium abundance employing the semi-analytical method introduced in Jlfjj . 
We remark that one can use numerical codes [11] for a proper treatment of BBN and accurate computation 
of the light nuclei abundances. However, the final density of He is very weakly sensitive to the whole 
nuclear network yfl. Therefore, in the present investigation we shall employ the semi-analytical treatment 



of H, computing the Helium abundance to a very good approximation avoiding sophisticated computer 
softwares. See also I12I1 for a recent example of a published work in which the same semi-analytical 
method was used to constrain the higher dimensional Planck mass in a brane model. Our results show 
that BBN requires the models considered in the present work to be only slightly different from the usual 
Einstein's general relativity, whereas dark matter consideration alone does not seem to constrain this class 
of new gravity theories due to the degeneracy in paramater space of the underlying particle physics models. 
However, from the BBN consideration we give a precise range for the allowed values of the power n and 



confirm the result of II13H using a different approach based on physics of the early universe. Finally, we 
remark at this point that according to our findings, certain scenarios that require a value of n considerably 
different than one cannot work. Furthermore, the models that satisfy our constraints do not lead to the late 
cosmic acceleration. 

Our work is organized as follows. The article consists of five sections, of which this introduction is 
the first. The modified gravitational model is described in the next section. The analysis based on BBN is 
discussed in section 3, while the investigation based on WIMP dark matter is presented in section 4. Finally 
we conclude in the last section. 



II. THE MODIFIED GRAVITATIONAL MODEL 

Here we shall present the model of f{R) gravity that will be discussed in this paper, and we shall 
summarize the basic formulas following [8]. The model is described by the action 

S= 2^J d V=S/(*) + S <" [#Mv, <M > (1) 
where G is Newton's constant, K 2 = &71G, and S m is the action of the matter field, <p m . Varying this action 
with respect to the metric we obtain the field equations for gravity, which generalize the usual Einstein's 
equations, 

f'Ruv - \fgnv ~ VpVvf + gpvOf = K% v , (2) 

where T^ v is the energy-momentum tensor for the matter, and a prime denotes differentiation with respect 
to R. For the gravity part we consider the spatially flat Robertson-Walker (RW) line element 

ds 2 =dt 2 -a(t) 2 (dx 2 + dy 2 + dz 2 ), (3) 

while for the matter part we consider a cosmological fluid characterized by a time-dependent energy density 
p(t) and pressure p{t) 

T$ = diag(p,-p, -p, -p). (4) 
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The — component of Q gives 



a . 1 d ,. . 9 
3-f'--f + 3-f"R = K 2 Pl (5) 
a 2 a 



while the £ — i components give 

'*+d^)f + lf- 2 h -f"R - f'R 2 - f'R = K 2 p. (6) 
a a J 2 a 

Here a dot denotes differentiation with respect to the cosmic time t. In addition to the above cosmological 
equations, we have as usual the energy conservation law 

p + 3- a {p+p)=0. (7) 

We restrict ourselves to the models of the form 

/(*)=(£)". <s> 

where A is a constant, A ~ M 2 p 2 ^", with M p = 1.22 x 10 19 GeV the Planck mass. The power n is the unique 
parameter of this class of modified gravitational models, and n = 1 corresponds to the usual Einstein's 
theory. Since we are interested in the physics of the early universe, we consider the radiation dominated 
era in which p = p/3, and p ~ a~ 4 . Searching for a power law solution for the scale factor, a(t) ~ t a , the 
cosmological equations determine the unknown power a in terms of n as follows: 

« = (9) 

Notice that when n = 1 we recover the known result a(t) ~ t l l 2 for the radiation era in the usual Einstein's 
general relativity. Then using the expression for the energy density 

P = ^8*T\ (10) 



one obtains the relation between time and temperature 

T _( 15 y/ 4 1/4 uf 

\4K*gJ 8a ?«A«/ 2 ' 
where counts the relativistic degrees of freedom for energy density, and 



(11) 



a =6 za a Ia — — - . (12) 



2(1 -2a) 1 - 2 " ' 

Since the quantity g a must be positive, the allowed range of a is 0. 155 < a < 1/2, and a = 1/2 corresponds 
to Einstein's theory. Finally the Hubble parameter is given by 



i 



i 



which generalizes the usual formula H{T) ~ T 2 , valid in the standard cosmology of Einstein's general 
relativity. 
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III. BIG-BANG NUCLEOSYNTHESIS 

In this section we briefly review the sequence of basic events leading to the synthesis of primordial He- 
lium during the early stages of the expansion of the universe in the standard cosmology based on Einstein's 
general relativity, following Jioll . Then we shall present the corresponding discussion and our results for 
the modified gravity case. Since all the parameters are fixed and the power n is the only free parameter of 
the model, our discussion will allows us to determine the allowed range for n. 

When the rates of the weak interactions keeping baryons in chemical equilibrium with leptons become 
comparable to the Hubble parameter, the neutron fraction X = n n /(n n +n p ) is frozen at some value X(T ~ 0) 
to be determined below, where n n and n p are neutron and proton number density respectively. Once the 
temperature has fallen below about 1/25 of the Deuterium binding energy, the Deuterium bottleneck opens 
up, and nearly all of the original neutrons present at the decoupling time are captured in A He. Taking into 
account the neutron decay, the final helium mass fraction is given by 

Y 4 ~ 2exp(-t c /z)X(T ~ 0) (14) 

where % = 885.7 ±0.8 sec [14] is neutron's lifetime, and t c ~ 3min is the capture time at which neutrons 



are captured into Deuterium. 

Now we discuss how to compute X(T ~ 0) and t c . To this end, we employ a semi-analytical method (see 
e.g. Iioll . 115]) which is sufficiently accurate and very useful, since the physics is very transparent and the 
dependence of the abundances on input parameters can be clearly worked out without sophisticated com- 
puter softwares. Indeed, the abundance of 4 He is very weakly sensitive to the whole nuclear network oO, and 
therefore a sufficiently accurate result can be obtained without using computer codes. Our semi-analytical 
approach relies on the work of [ 10], and we shall not present here in detail all relevant formulas, as they are 
quite involved. To compute X(T ~ 0) we need to integrate the following rate equation 

^ = X pn (t){\ -X(t)) - KpXit). (15) 

Here we denote by X pn the rate for the weak processes to convert protons into neutrons and by X np the rate 
for the reverse processes that convert neutrons into protons. These rates are time dependent because of their 
temperature dependence. The rate X np is the sum of the rates of three processes 

Xn P = X(v+n^p + e-)+X(e + + n^p + v)+X(n^p + V + e-), (16) 

each of which is computed using standard field-theoretic techniques. After a few simplifications the rate 



X np is computed as follows jl0(1 



V v) - ( ^T) (12 + 6J+J 2 ), (17) 
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where y = Am/T with Am = m n — m p = 1.29 MeV being the neutron-proton mass difference. The detailed 
balance relation gives 

h P n(y) = e- y Kp(y)- (18) 

We now rewrite (fT5l) in terms of y instead of time as 

= T {Xp " m - x (y))-^ P (y)x(y)) , (19) 

where dt/dy can be computed using the definition of y, y = Am/T, and the fact that T/T = —H. With the 
initial condition X(y = 0) = 1/2, the rate equation for X(y) can be integrated numerically, and from the 
graphical solution one can compute X(T ~ 0), which we denote by X. 

Finally, let us add a few words regarding the capture time. The bottleneck opens up when the main 
reaction converting Deuterium into heavier elements 

D + D^T + p (20) 

become efficient. First we introduce 

z= £ f, (2D 

where Eq = m p + m„ — m D = 2.23 MeV is the Deuterium binding energy. The another quantity that is 
important in estimating capture time is the Deuterium abundance, Xq = nj)/{n n +n p ), where rip is the 
Deuterium number density. From the Saha equation Xd is given by 

X D = 2.8 x \0- l4 7 lw T^ / e 2 v exp(z)X p X, (22) 

where X p = n p /(n n + n p ) is the proton fraction. In the above formula TMeV is the temperature in MeV units, 
and we have parametrized the baryon-to-photon ratio by 

t]io = 10 10 x-, (23) 

fly 



where we use the observational value nt/riy = 6. 1 x 10 10 from WMAP [18]. The condition that determines 



the temperature (or time) at which the Deuterium bottleneck opens up reads as follows(for more details we 

n 

refer the reader to [10]). 

2X D R DD ~1, (24) 

where 

R DD = ^L <G V> n h = 2.9 x 10V 4/3 ec/>(-1.44z 1/3 ) (25) 
dz 
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and < ov > is the thermal average of the relevant cross section times relative velocity, which is a function 
of z. Then one can solve (1241) with respect z to get t c , and then from (fT4l one can finally obtain Y4. 

We now consider the constraints on f(R) gravity coming from BBN by computing the Helium mass 
fraction at the conclusion of the BBN. To this end, we will follow H for standard cosmology as described 
above, and the modifications will be done by adopting the relation (fTTT) instead of the standard one. 

In the modified gravity model, the basic physics governing the details of primordial nucleosynthesis 
remains the same, and the only thing that is different now is the new time-temperature relation (fTTT >. from 
which one obtains 
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47T 3 10.75 



l/4a M l/2a 



A 1 2 \Am) 



l/a 



l/4a 



(— 



l/4a M l/2a 



l/a 



l/4a 
'a 



(26) 
(27) 



.37/ AV2 

In (|261)-(|271) we have taken into account the appropriate value for at the relevant temperature. Finally, the 
4 He mass fraction is still given by (fT4l . but now both freeze-out abundance X(T ~ 0) = X and capture time 
t c are modified due to the new time-temperature relation. 



8 = l-n 


X 


f t (sec) 


*4 





0.1529 


176.76 


0.2504 


1(T 5 


0.1528 


176.94 


0.2503 


IO- 45 


0.1526 


177.34 


0.2499 


10- 4 


0.1521 


178.59 


0.2486 


IO- 39 


0.1519 


179.06 


0.2482 


10 -3.8 


0.1516 


179.66 


0.2476 


io- 3 - 7 


0.1513 


180.42 


0.2468 


10 -3.6 


0.1509 


181.37 


0.2459 


1Q -3.5 


0.1504 


182.58 


0.2448 


IO" 34 


0.1497 


184.12 


0.2433 


IO" 33 


0.1489 


186.07 


0.2414 


1Q -3.2 


0.1479 


188.54 


0.2391 


IO" 31 


0.1466 


191.71 


0.2362 


IO" 3 


0.1450 


195.77 


0.2326 


io- 2 - 9 


0.1430 


201.01 


0.2280 



Table I: Helium abundance, capture time, and freeze-out neutron mass fraction for several values of 8 = 1 — n. 



First we integrate (fl9l ) with the initial condition X(y = 0) = 1/2 to obtain X (it is enough to evaluate X at 
y = 15 because it freezes out). Then we use the condition (l24l) to compute t c . Note that now the function t(T) 
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is the one predicted by the new gravitational model. Finally we compute Y4 from (fT4b for several values of 
a (or n). Our results are shown in Tableland Fig. 1. The results show that F 4 is very sensitive to 1 — n. The 
freeze-out abundance decreases, while the capture time increases, as 1 —n increases (or a decreases). These 
both effects bring about decrease in Y4 as 1 — n increases. The reason is as follows. The decrease of the 
freeze-out abundance is due to the decrease of the freeze-out temperature, which is calculated by equating 
the Hubble parameter (PT31) with the weak interaction rate F = X np {T) + X pn {T) ~ r 5 jl6ll . The increase of 
the capture time is mainly due to the temperature-time relation: For fixed temperature, the corresponding 
time increases as a decreases. 

We have interpolated the computed helium abundances as a function of n (Fig. 1) and obtained the lower 
bound of n in comparison with the observational data. The allowed range for n is 

l-n< 0.00016 or equivalently a > 0.49992 (28) 

according to the observational constraint jl7ll 

Y 4 = 0.2516 ± 0.0040. (29) 

We see that according to our results the power n, compared to the value found in [8], should be pushed even 
closer to unity. This implies that the predicted baryon asymmetry within the framework of gravitational 
baryogenesis using this class of models becomes now too low. Furthermore, several existing scenarios that 
require a value of n considerably different than unity cannot work. 

IV. WIMP DARK MATTER 

Recent cosmological observations Jisll have established the allowed range of the normalized density of 
cold dark matter in the universe 

0.075 <a cdm h 2 < 0.126. (30) 

In the present section we assume that the role of cold dark matter in the universe is played by weakly- 
interacting massive particles (WIMPs) proposed by physics beyond the standard model. For concreteness 
one can think of the lightest neutralino. The discussion to follow is a model-independent one, in which we 
have considered a generic WIMP assuming that its mass is ( 100 — 500) GeV, and that its typical cross section 
of the relevant processes in which it participates is not very different than o ~ aj m /M^ w , where a em = 1/137 
is the electromagnetic fine structure constant, and M ew ~ 100 GeV is the electroweak scale. The relic density 
of the dark matter particle depends on its mass m, its annihilation cross section Go, and finally on the power 
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n characterizing the gravitational model. Since BBN has already determined the allowed range for n, we fix 
it to a given value and therefore the cold dark matter constraint Q.h 2 ~ 0.1 gives a certain relation between 
the WIMP mass and the annihilation cross section. 

The evolution of the number densit y n of the dark matter particle in an expanding universe is determined 



by solving the Boltzmann equation 



isit y n o 

fly 



n + 3Hn = - < OV > (n 2 -n 2 EQ ), (31) 



n E Q=l 7 ^f(E). (33) 



where H is the Hubble parameter, heq is the number density at equilibrium, v is the relative velocity, and a 
is the total annihilation cross section. The thermal average of the total annihilation cross section times the 
relative velocity < av > is given by 

< ™>=^Jm^ fiEl)nEl)av < (32) 

where f(E) is the fermion distribution function, f{E) = 1/(1 + exp(E /T)). Finally the number density at 
equilibrium is given by 

d 3 p 

Now it is convenient to introduce new variables, namely dimensionless quantities 

m 

x = -, (34) 
n 

Y = -, (35) 
s 

where T is the temperature and s is the entropy density 

2?r 2 

S = K ^ T3 (36) 
with h* being the number of relativistic degrees of freedom for entropy density. Assuming entropy conser- 
vation, the Boltzmann equation can be written down equivalently as follows 

§--£<«» XT*-*). 07) 

The yield at equilibrium Yeq for non-relativistic (cold, i>3) relics is given by the approximate expression 

45 fTi\ l l 2 x 3 / 2 exp(-x) 

where g = 2 is the spin polarizations of the dark matter particle. In standard cosmology during the radiation 
dominated era, the Hubble parameter as a function of the temperature is given by H(T) = \ .61g l J 2 T 2 /M p . 
Parameterizing < av > as 

< ov >= Oqx~ 1 (39) 
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the Boltzmann equation takes the final compact form 



^ = -A*-'- 2 (F 2 -Fj e ), (40) 

where A is a constant given by 

k = ( X< ° V > S ) = 0.264(h/gl /2 )M p mo . (41) 

V H { m ) J x=l 

We can obtain an approximate analytical solution of Boltzmann equation by the following arguments. Ini- 
tially, for large temperatures the annihilation rate is larger than the expansion rate of the universe and the 
WIMP abundance follows the equilibrium abundance. At some point x/ the annihilation rate becomes com- 
parable to the expansion rate and the dark matter particle decouples from the thermal bath. For x 3> x/ we 
can neglect the Yeq term in the Boltzmann equation. Then the equation can be easily integrated and the 
solution Y oa = Y(x = °°) is given by 

Y 00 = l -±±x l ;\ (42) 



where the freeze-out temperature Xf is given j 1611 by 



x f = ln[0.038(l + l)(g/gl /2 )M p mo } (43) 
I + ln{ln [0.038 (/ + l)(g/gl /2 )M p m(Jo}}. 

After having integrated the Boltzmann equation for Y(x), then the relic abundance for the dark matter 
particle is given by 

a cdm h- = ^^, (44) 

Per 

where Tq is the today's temperature. Here we make use of the following values: 

T = 2.73K = 2.35 x 10~ 13 GeV (45) 
h*(T ) = 3.91 (46) 
p cr /h 2 = 8.1 x 10~ 47 GeV 4 (47) 

So far we discussed the case of the standard cosmology. Now let us take into account the modification 
of gravity. Taking into account a ~ t a and time-temperature relation (fTTT> . one can express the Hubble 
parameter as 

H = H a (m)x-^ a , (48) 
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where 



i 



aA2 f 4k g 



*4a ]\/f2a \ 15 



r3„ \ 



1 



H a {fn) = — r ( I . (49) 



For a = 1/2, this reduces to the usual H(m) parameter for the standard cosmology llql 



Ha=i/iH = l.67gl /2 m 2 /M p =H(m). (50) 
For x > 3 the temperature dependence of the annihilation cross section is parameterized as 

< av >= o x~ l , (51) 

where / = corresponds to s-wave annihilation, / = 1 to p-wave annihilation, etc. Then the Boltzmann 
equation for the abundance of dark matter becomes 

^ = -lx- ! -HY 2 -Yl Q ), (52) 

where 

- = ( = PfirX = 0.264(K/g^)M pm o . (53) 
V H a (m) J x=l H a {m) 

Here we introduced new parameter I and do to clearly show the effect of the modification of gravity on 

Boltzmann equation: 

/ = / + (2-i), (54) 

H(m) 
H a (m) 

It follows that I = I and do = Cb for a = 1/2. It is clear that the modification coming from f(R) is 
implemented entirely by the correction on these two parameters. By comparison to (|401)-(|4TI). one can 
easily see that the Boltzmann equation (1521 together with (1531 exactly corresponds to one of the standard 
cosmology with the averaged product of annihilation cross section and velocity 

< av >= &qx~ 1 . (56) 

Since the Boltzmann equation has exactly the form as in standard cosmology, one would get the same results 
as in the standard case, but with the replacement, Go — > Go and / — > Iq. 

The most important quantity in estimating the relic density is x/, which is the time when Y ceases to 
track Yeq, or equivalently, when Y — Yeq becomes of order Yeq. This quantity is computed by (1431 as 

x f = ln[0.038(/ + l)U/^ /2 )M p ma ] (57) 
f+ ^) ln{ln[0.038(f+ l){g/gl /2 )M p mo ]}. 
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Then the present yield 



Yoo and relic density Q. c d m h 2 are given by 



3.79(1+ 1)/ 



(58) 



1 /2 

[K/gJ )M p m&o 



Q. cdm h 2 = 1.07 x 10 9 



(f+l)x' f +1 GeV 



(59) 



{h,/gl /2 )M p d 



Since the cosmic temperature during the period of interest is T ~ (a few)GeV, we set h* = g* ~ 100. Eq. 
(TSTT ) - ( T59b are one of the main results of this section. It is remarkable that they are expressed as an analytical 
function of a through /, Go and Xf. Thus once I, Go and m together with a are given, the relic density is 
directly obtained from (l59l . We will use s-wave approximation (/ = 0), so I = 2 — I/a. In this case / is 
negative because a < 1/2, but since we know from the consideration in the previous section that BBN 
allows tiny deviation of n from 1 in ([8]), we assume that / + 1 > in (l57l)-(|59"l). 

We are finally in a position to present our numerical results in figures, showing the relation between the 
annihilation cross section and WIMP mass. The dark matter abundance is a function of three parameters, 
namely n,m, Go- If the power n is fixed according to the BBN results, and we impose the cold dark matter 
constraint 0.075 < D. cc i m h 2 < 0.126, it is possible to obtain a certain relation between the WIMP mass m 
and its annihilation cross section Go. Our results can be shown in the figures 2 and 3 below. In particular, 
in Fig. 2 we show the annihilation cross section as a function of the WIMP mass for fixed values of n 
(corresponding either to general relativity or to the new gravitational model for the range determined from 
BBN), and for the upper limit Q. c d m h 2 = 0.126. In Fig. 3 we show Go as a function of m for fixed values of 
n and for the lower limit Q. C( i m h 2 = 0.075. From these figures one can see what the lower bound (Fig. 2), 
and upper bound (Fig. 3) of the annihilation cross section should be for a given WIMP mass. 



In the present work we have studied primordial Big-Bang nucleosynthesis and WIMP dark matter in 
a class of modified gravitational theories. This class of gravitational models predict a novel expansion 
law for the early universe. For BBN we have employed a semi-analytical computation in which the basic 
physics is quite transparent. For WIMP dark matter we have given a model independent discussion applying 
the usual treatment found in standard textbooks or reviews. Concerning BBN, by comparing the theoretical 
predictions to the available observational data we were able to put bounds on the unique parameter appearing 
in this class of modified gravitational theories. We have found that the models considered in the present 
work are allowed to be only slightly different from the usual Einstein's general relativity. In the dark matter 



V. CONCLUSIONS 
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section we have obtained an analytical expression for the cold dark matter abundance as a function of 
n. After that we fixed the power n according to the BBN results, and we have shown in figures how the 
annihilation cross section and the WIMP mass are related in order that the cold dark matter constraint is 
satisfied. The predicted baryon asymmetry within the framework of gravitational baryogenesis using this 
class of models [a] is too low for our BBN range of n, and therefore this mechanism for baryon asymmetry 
does not seem to be consistent with BBN constraints. Finally we remark in passing that the models that 
satisfy our bounds do not lead to the late cosmic acceleration. 
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Figure 2: Annihilation cross section versus WIMP mass for CDM abundance Q. cc i m h 2 = 0.126. Shown are n = 1 
(solid), n = 1 - 1(T 4 (dashed), and n = 1 - 2 x 1(T 4 (dotted). 



